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related to non-degeneracy of heavy Majorana neutrinos? 
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Abstract 

We propose a scenario that the mass splitting between the first generation of the heavy Majorana 
neutrino and the other two generations of degenerate heavy neutrinos in the seesaw framework is 
responsible for the deviation of the solar mixing angle from the maximal mixing, while keeping the 
maximal mixing between the tau and muon neutrinos as it is. On top of the scenario, we show that 
the tiny breaking of the degeneracy of the two heavy Majorana neutrinos leads to the non-zero 
small mixing angle Ues in the PMNS matrix and the little deviation of the atmospheric neutrino 
mixing angle from the maximal mixing. 
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Thanks to enormous progress in solar, atmospheric and terrestrial neutrino experiments, 
we have now the robust evidence for the existence of neutrino oscillation which provides 
a window to physics beyond the standard model (SM). Until now, while the atmospheric 
neutrino deficit still points toward a maximal mixing between the tau and muon neutrinos, 
however, the solar neutrino problem favors a not-so-maximal mixing between the electron 
and muon neutrinos. There have been many attempts to explain the origin of the deviation 
of the solar mixing angle from the maximal mixing. Surprisingly, it has recently been noted 
that the solar neutrino mixing angle 6soi required for a solution of the solar neutrino problem 
and the Cabibbo angle 9c reveal a striking relation 9soi + — f , which is satisfied by 
the experimental results within a few percent accuracy, 6soi + dc = 45.4° ±1.7° . This 

quark-lepton complementarity (QLC) relation has been simply interpreted as an evidence 
for certain quarU-.epton symmetry or quarU-.epton uuiflcation ^ shown in Refs. 
But, it can be an accidental phenomenon as pointed out in Ref. [6|,13|. Thus, it is worthwhile 
to find the possible alternatives to the grand unification origin of the deviation of the solar 
mixing from the maximal mixing. 

In this letter, we propose a scenario that the mass splitting between the first generation of 
the heavy Majorana neutrino and the other two generations of degenerate heavy neutrinos 
in the seesaw framework is responsible for the deviation of the solar mixing angle from the 
maximal mixing, while keeping the maximal mixing between the tau and muon neutrinos 
as it is. The maximal atmospheric neutrino mixing and the smallness of f/es may be the 
trace of the original "bi-maximal" mixing which is presumably supposed to be achieved by 
some underlying flavor symmetries, and thus the best possible approach to the problem is to 
start in the limit of the maximal mixing with [/es = 0, and understand how the deviation of 
the solar mixing from the maximal is realized. In our scenario, the primitive "bi-maximal" 
neutrino mixing is generated only from the neutrino Dirac Yukawa matrix by taking a 
diagonal form of three degenerate heavy Majorana neutrinos in a basis where the charged 
lepton mass matrix is real and diagonal. As will be shown, the deviation of the solar mixing 
can then be generated from breakdown of the degeneracy of the heavy Majorana neutrino 
masses between the first and the other two generations. The main point in this scenario 
is that the deviation can be expressed in terms of the ratio between two heavy Majorana 
neutrino masses. On top of the scenario, we will also show that the tiny breaking of the 
degeneracy of the two heavy Majorana neutrinos will lead to the small mixing angle ^13 in 
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the PMNS matrix and the very small deviation of the atmospheric neutrino mixing angle 
from the maximal mixing. 

Before proceeding to our scenario, we wish to motivate one scheme that leads to exact 
"bimaximal" mixing in the framework of the seesaw mechanism. We study in a basis where 
the charged lepton mass matrix is real and diagonal. The light neutrino mass matrix My 
diagonalized by f/bimax is given through the seesaw mechanism by 

M, = Ml M^' Mr), 

= f^bimax M^"^^ f^bimax i (l) 

where Mr, = Yr) v/ \/2 with electroweak vacuum expectation value v and the neutrino Dirac 
Yukawa matrix Yd, and Mn is a mass matrix of heavy Majorana neutrinos. The mixing 
matrix f/bimax denotes the "bi- maximal" mixing matrix 



f/bimax — U2-i (y^) Ui2 
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Then, the "bi-maximal" mixing can be achieved by one of the three possible ways as follows: 

• Taking Yd diagonalized by f/bimax and M^ = M ■ I with the identity matrix / and a 
common mass scale M. 

• Taking Yd = y ■ I and Mr diagonalized by f/bimax- 

• Taking "bi-maximal" mixing pattern from the combination of the nontrivial Yd and 

For the third case, there may exist various origins of the deviation of the solar mixing 
depending on possible combinations, and some of which have been discussed before [9j . For 
the other two cases, the modification of the trivial sectors proportional to the unit matrix 
can be in charge for the origin of the deviation from the maximal mixing. However, since the 
second case may lead to the undesirable deviation of the atmospheric mixing aside from the 
deviation of the solar mixing as one can easily see, we only focus on the first case in this letter. 
In the first case, the "bi-maximal" mixing can be achieved by taking the symmetric matrix 



Yd with specific form. As an example, we present a detailed model of Yd leading to the "bi- 
maximal" mixing, while keeping Mn = M ■ I based on the discrete symmetry A4 (g) Z2 Q| . 
Let the three families of leptons and singlet heavy neutrinos be denoted by (z/j, li)^, liji, Nipi 
for i = 1,2, 3. In this convention, luljR and VnNjfi are Dirac mass terms for charged leptons 
and neutrinos. Under the discrete symmetry A^^® Z2, the 3 families of leptons transform 
as {Pi, li)L ~ (3, +), NiR ~ (3, +), liR ~ (1, -), (1', -), (1", -). We introduce Higgs scalar 
sectors consisted of seven Higgs doublets $i ~ (1, — ), (1', — ), (1", — ), (p ~ (1,+), cr^ ~ 
(3, +). From the assignment, the ® Z2 invariant Dirac Yukawa interactions for charged 
lepton sector, lijliR^j, leads to a diagonal mass matrix with 3 independent entries each as 
shown in Ref. For the mass matrix of the heavy Majorana neutrinos, we can take 

MNiFtNiR with common mass scale M because of ^4 symmetry, i.e. 3 x 3 ~ 1. The Dirac 
Yukawa matrix for the neutrino sector, which is invariant under A4 (S> Z2 and diagonalized 
by the "bi-maximal" mixing matrix, can be obtained from the interaction Lagrangian as 
follows: 

Yd = hiuiNi + iy2N2 + iy3N3)(l) 

+ h2{iyiN2cr3 + iy2N3cri + U3Nia2) 

+ h3{NiV2<y-i + N2V3(7i + N3Via2) + h.c. (3) 

In order to achieve the symmetric form of the Dirac Yukawa matrix, we require h2 = h^. The 
vacuum expectation values for the neutral components of Higgs sector cr° can be determined 
by the Higgs potential invariant under A/^, 
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V = m ajai + -\i{aja^ 



+ A2(cr}o-i + iu o\a2 + ujalp':^{o\ax + uja\o2 + 0-30-3) 
+ A3[(cTjo-3)(cT;^a2) + {p\oi){a\o-3) + {p\a2){(j\ox)\ 

+ {^A4[(ata3)^ + {a),a,f + {o\a2f\ + /..c.} , (4) 



as 



where u; = e^'^l'^ . Taking < o\ >= and < 0-° >=< 0-3 >= v with v — y 2A1+A2+A3+A4 
well as non-vanishing < cfp > for the Higgs sector 0, we can achieve the final form of the 
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Dirac Yukawa matrix given as follows, 



abb 
b a 
b a 



(5) 



Defining Y^'^^ = diag{x, y, z), the neutrino Dirac Yukawa matrix Yd diagonalized by C/bimax 
is generally given in the symmetric matrix form by 



V — TT V^^^S TjT 

J^D — (^bimax d <^bimax ■ 



(6) 



Here, we consider the case of nonzero values for x and y, which is crucial to our purpose. 

In order to achieve the observed deviation of the solar neutrino mixing from the maximal 
mixing, we take into account the mass splitting between the first generation of the heavy 
Majorana neutrino and the other two degenerate ones, for which the mass matrix is given by 
Mr — M^^^ — (Ml, M2, M2), which results from the breaking of A4 in the heavy neutrino 
sector and refiects separation of Nir ~ Nir{1) -/V(2,3)r(2) under S3 symmetry. Then, the 
light neutrino mass matrix Mj, is presented as follows: 
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= ^bimax-^i/t^jimax i 

where the mass matrix M'^ is given by 

/ 



Mo 



12 



X 



\ 



bimax 



(7) 



M' = 



(Ml + M2) 



2M1M2 



(Ml - M2) 



2M1M2 

2&iM^-M,) ^(Mi + M2) 

ij 



(8) 



Then, the matrix M^ can be diagonalized by UuiO), and after diagonahzing M^, we can 



5 



obtain the mixing angle 9 and three neutrino mass eigenvalues as follows: 

^ [(c'x' + sV)(Mi + M2) + 2csxy{Mi - M2)], 



2M1M2 



^ [(sV + cV)(Mi + M2) -2csxy(Mi-M2)], (10) 



2M1M2 

2 



" M2 ' 

where c = cos^, s — sinO. Comparing the mixing matrix Ui2{9) with C/i2(7r/4) in Uumax, 
we can get the solar mixing angle Ogoi which deviates as much as the value of 9 from the 
maximal mixing. Note that the value of 9 should be negative in order to achieve the desirable 
deviation of the solar neutrino mixing. We can argue that the generation of the mixing angle 
9 due to the splitting between Mi and M2 in seesaw mechanism may be the origin of the 
deviation of the solar mixing angle from the maximal mixing in the case of nonzero x and y. 
However, since we do not have yet any information on the values of Mi and M2, we cannot 
immediately test whether the difference between Mi and M2 is really compatible with the 
deviation of the solar mixing angle from the maximal mixing, but we can make numerical 
estimate for the size of the ratio of Mi to M2, which accommodates the deviation of the solar 
mixing based on the experimental results for the neutrino oscillation. From the numerical 
results, we can also predict the magnitude of the effective Majorana neutrino mass niee, 
which is the neutrino-exchange amplitude for the neutrinoless double beta decay. 
For our purpose, let us define two parameters k and cu as follows: 

^=x^ "=M2- ^''^ 
Then, the expressions for 9 and m^,. are given as follows, 

tan 2^ = 7 — (12 

(l-«2)(i+a;)' ^ ^ 

''2, „2,2 



[(c^ + sV^)(l + a;) + 2csk{uj - 1)], 
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2 
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In addition, the effective Majorana neutrino mass viee is presented by 

mee-^[{l + ^?+U;{l-K)\ (14) 



As shown in Eq. fll2|) . the non- vanishing value of the mixing angle 9 can arise when uo is 
deviated from one, which indicates the splitting between Mi and M2. In fact, the present 
experimental results are not enough to determine all the parameters introduced. But, if we 
fix one neutrino mass eigenvalue by hand, we can determine several independent parameters 
as well as the magnitude of from Eqs. (|12ll3ll4j) . For our numerical calculation, we set 
the parameter 9, Aml^ and Aml2 to be 13°, 8 x 10~^ eV^, 2.5 x 10~^ eV^, respectively. 
Those numbers correspond to the best fit values for the measurements of the deviation of 
the solar mixing angle from the maximal mixing, the mass-squared differences of the solar 
and atmospheric neutrino oscillations, respectively. By fixing m,^^ as an input parameter, 

2 2 

we can determine the parameter set (fi:,a;, -^) for normal hierarchy m^-^ < m^^ < m^^ 
through the relations (jllll2ll3p . 

TABLE I: All numbers corresponding to the mass parameters are given in the unit eV for normal 
hierarchy. 

m^l (input) K UJ ^ ^ TUee 

0.005 1.298 0.772 0.003 0.051 0.009 

0.01 1.118 0.897 0.006 0.052 0.013 

0.05 1.006 0.994 0.025 0.071 0.051 

0.1 1.002 0.998 0.050 0.112 0.101 



In TABLE I, we present our numerical results for normal hierarchy. From the TABLE I, 
we can see that the values of k and u approach to one as m,^^ increases up to of order 0.1 eV, 
and one needs fine-tuning to obtain the parameter set satisfying the relations above for the 
case of such a large m^-^ ~ 0.1 eV. As mj,^ goes down, the value of k, rapidly increases whereas 
that of uj decreases. We can also predict the size of the amplitude of the neutrinoless double 
beta decay function of rriy-^ , which is presented in the last column of TABLE I. If the 

neutrinoless double beta decay wih be measured in near future, we wih be able to determine 
three neutrino mass eigenvalues and the parameters introduced in Eqs. p2ll3ll4j) . For 
inverted hierarchy iriy^ < m^^ < m^,^, the numerical results are presented in TABLE II. In 
this case, contrary to the normal hierarchical case, we take m^^g as an input. 

Next, to generate non- vanishing f/es, on top of the above scenario, we consider an in- 
teresting possibility that the breaking of the degeneracy between the second and the third 
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TABLE II: 
hierarchy. 



All numbers corresponding to the mass parameters are given in the unit eV for inverted 



1,3 (input) 


^ Ml M2 "^ee 


0.005 


1.672 0.585 0.011 0.010 0.041 


0.01 


1.569 0.630 0.013 0.013 0.042 


0.05 


1.137 0.881 0.029 0.051 0.066 


0.1 


1.044 0.959 0.052 0.100 0.109 



generation masses in the heavy Majorana neutrino sector, i.e., Mr — diag{Mi, M2, M3), can 
be an origin of the generation of non- vanishing C/gs. We remark that the value of C/gs goes 
to zero in the limit of M2 — M3 in this scenario. The effective light Majorana neutrino mass 
matrix is given by 
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Assuming that the mass splitting between M2 and M3 is small enough to accommodate 
the tiny f/e3, the mixing matrix, which diagonalizes the neutrino mass matrix M^,, can be 
approximately given by 
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(17) 



where the mixing angle a corresponds to the (1,2) rotation of 2 x 2 submatrix of M'^. The 
mixing angle a is presented by 

2K{l-uj-e) 



tan (J ~ 



(18) 



(l-fi:2)(l + cu + £)^ 

where e — M1/M3, and u!,k are given earlier. This mixing angle a is responsible for the 
deviation of the solar mixing angle from the maximal mixing. We note that non- vanishing 
value of a is possible even when u; — 1, i.e. {Mi — M2), but this case is undesirable because 
it leads to negative a which positively contributes to 9i2. The mixing angle a is zero when 



uj + e = 1, but it corresponds to the large hierarchy among three heavy Majorana masses, 
which is far beyond our purpose. The mixing elements S and r] are given by 



1 1 



where ci and C2 are presented in terms of three light neutrino mass eigenvalues and the 
parameters K,u,e. Then, the mixing element Ues and the deviation of the atmospheric 
mixing from the maximal mixing are simply presented in terms of cr and f] as follows, 

\Ue3\ ^ l\S-v\, (20) 
5 sin ^23 - ^{S + V)- (21) 
Imposing the bound on {Ue^l of CHOOZ experiment, If/gsl < 0.2, and the result of sin^ ^ 
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from atmospheric neutrino data, sin^ ^23 = 0.44(1^022) at 2a [1^, we can determine the 
allowed regions of the ratio M2/M3. In TABLE III, we present the numerical results for 
the ratio M2/M3 and the prediction for the bound on {Uesl- The second and third columns 
correspond to the normal hierarchical case, whereas the fourth and fifth columns to the 
inverted hierarchy. We find that the result for 5sin^^23 constrains M2/M3 more severely 
than the bound on \Ue3\ for m,^^^^^ < 0.05 eV. But for m^^^^^ ~ 0.1 eV, both 5sin^6'23 and 
If/esI from neutrino data severely constrain the allowed region of M2/M3. The values in the 
columns for {Uesl indicate the predictions for the upper bound. As shown in Table III, the 
allowed region for M2/M3 gets narrowed as ^i/^g) increases, and it becomes nearly one for 
^vi(3) > 0.1 eV. This implies that such large values of m^^;^^ lead to moderately degenerate 
light neutrino spectrum realized by almost degenerate heavy Majorana neutrinos. 

Finally we note that there could be radiative corrections to neutrino mass matrix which 
can lead to some modification of our results. However, non-negligible renormalization effects 
can be expected only in the case of degenerate light neutrino spectrum. The numerical results 
for m;^j(3) = 0.1 eV in the tables may be significantly modified due to possible renormalization 
effects, but the detailed investigation on the renormalization effects is not our main interest 
in this work and we will leave it for the future work. 

In summary, we have proposed a scenario that the mass splitting between the first gen- 
eration of the heavy Majorana neutrino and the other two degenerate ones in the seesaw 



TABLE III: The numerical results for the ratio M2/M3 and the prediction for the bound on \Ue3 
for the normal hierarchical case and the inverted hierarchical case. 



mi.^(3j (input) 


M2/M3 \Ue3\ 


M2/M3 \Ue3\ 


0.005 


0.55 ~ 1.29 (< 0.015) 


0.44 ~ 1.36 (< 0.025) 


0.01 


0.72 ~ 1.16 (< 0.0007) 


0.59 ~ 1.27 (< 0.015) 


0.05 


0.94 ~ 1.04 (< 0.086) 


0.93 ~ 1.04 (< 0.083) 


0.1 


0.98 ~ 1.01 


0.98 - 1.01 



framework is responsible for the deviation of the solar mixing angle from the maximal mix- 
ing, while keeping the maximal mixing between the tau and muon neutrinos as it is. Our 
scheme is based on the assumption that nature presumably started with "bi-maximal" neu- 
trino mixing and then it has been deviated somehow. We have considered the case that the 
"bi-maximal" mixing is achieved only from the neutrino Dirac Yukawa matrix by taking a 
diagonal form of three degenerate heavy Majorana neutrinos in a basis where the charged 
lepton mass matrix is real and diagonal. Allowing the mass splitting between the first and 
the other two generations of the heavy Majorana neutrinos, we could obtain the deviation 
of the solar mixing angle from the maximal. In addition, we have also shown that the tiny 
breaking of the degeneracy of the two heavy Majorana neutrinos leads to the small mixing 
angle 6^3 in the PMNS matrix and the very small deviation of the atmospheric neutrino 
mixing angle from the maximal mixing. 
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